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Abstract 

We report the OPE-based predictions for a number of lepton energy and hadronic mass 
moments in the inclusive semileptonic B^X c £u decays with a lower cut on the charged 
lepton energy. We rely on the direct OPE approach where no expansion in the charm 
mass is employed and the theoretical input is a limited set of underlying OPE parameters 
including m;, and m c . A Wilsonian treatment with a 'hard' cutoff is applied using running 
low-scale masses mg(/i) and kinetic expectation value This leaves for perturba- 

tive corrections only genuinely short-distance effects and makes them numerically small. 
Predictions are also given for the modified hadronic moments of the kinematic variable 
M x which is a combination of M\ and Ex- Measurement of such moments would allow 
a more reliable extraction of higher-order nonperturbative heavy quark parameters from 
experiment. 



*On leave of absence from St. Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg 188300, Russia 



1 Introduction and motivation 



The heavy quark expansion based on the local Operator Product Expansion (OPE) 
[1-3] allows to accurately calculate sufficiently inclusive decay probabilities, incorporat- 
ing bound-state and hadronization effects in terms of a limited number of physical heavy 
quark parameters (for details, see the review [4]). Inclusive B decay distributions - in 
particular those related to civ transitions - represent a portal to a precise determi- 
nation of these parameters [5,6], and lepton energy moments and moments of hadronic 
mass and/or energy are among the most interesting quantities. The main ingredients for 
their computations have been around for some time (see [2,7-13]). Several numerical 
applications can be found in the literature, for moments with or without cuts on the 
charged lepton energy (see [14-16] and references therein) often dictated by the reality of 
experimental measurements. The present study differs from and complements previous 
reports in a number of aspects. 

First, we make use of the robust approach advocated in [17, 18]; in particular, we 
do not invoke an expansion in l/m c which has plagued the reliability of many earlier 
applications of the heavy quark expansion. Moreover, in our approach only the heavy 
quark parameters relevant to inclusive decay rates are invoked. This reduces the number 
of new objects appearing at order 1/toq from six to two and eliminates poorly known 
non-local correlators. An extended set of experimental moments allows to constrain all 
relevant parameters provided model-independent bounds [19] are incorporated. 

We also rely on heavy quark parameters which are renormalized a la Wilson and 
depend explicitly on a 'hard' normalization scale \i which we set equal to 1 GeV. Primarily, 
this refers to the heavy quark masses m\,{}i) and m c (/i) and to the kinetic expectation value 
a4(a*)- On the theoretical side, this is necessary both to meaningfully assign them definite 
values and to apply exact heavy quark inequalities. On the practical side this renders 
perturbative corrections well-behaved and moderate in size. Absence of large higher order 
corrections is crucial for a meaningful extraction of the nonperturbative parameters. 

The direct use of the underlying set of heavy quark parameters reveals what a partic- 
ular moment is actually measuring or constraining, and simplifies considerably the task 
of estimating the theoretical accuracy, often the subject of controversial claims [20]. It is 
worth reminding [6] in this respect that the inclusive decays cannot depend on non-local 
correlators often appearing in other applications of the heavy quark expansion, and there 
is no way to constrain them directly studying only inclusive B decays. 

The important motivation behind the present study is to make available a code for 
evaluating the various distributions, which is not bound to specific technical assumptions, 
for instance the use of meson mass relations, and is flexible enough to allow a meaningful 
investigation of the theoretical uncertainty using different options (see Section 3). The 
detailed analytic expressions will be given in Ref. [21], while here we limit ourselves with 
presenting easy numerical recipes and a discussion of the main points. 

Finally, we present here predictions for the modified higher hadronic moments with 
a charged lepton energy cut, which can be measured in the experimental setup of the B 
factories. They will allow to pinpoint higher-order nonperturbative expectation values 
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with better accuracy and reliability. 



2 Theoretical setup 

All inclusive semileptonic B decay distributions are described by a few B decay structure 
functions Wi(qo; q 2 ), where g M is the four-momentum of the lepton pair; for massless leptons 
one has [8] 

d^dg 2 dg 32tt 4 V w *Ei) v ej \h ) 

{2q 2 Wl + [4E t (qo-E e )-q 2 ]w 2 + 2q 2 (2E e -q )w 3 } . (1) 

The OPE-based heavy quark expansion yields them in terms of (short-distance) quark 
masses and of the 5-meson expectation values of local heavy quark operators. Accounting 
for the latter through D — 6, the set of input parameters for our evaluation of the moments 
includes therefore m b , m c , /i 2 , /Iq, p% and p 3 LS . Hadronic mass moments kinematically 
depend also on the B meson mass M B for which the experimental value 5.279 GeV is 
employed. In practice, it appears in the combination A' = M B — m b . As mentioned above, 
the input heavy quark parameters depend on the normalization point \x for which we 
adopt 1 GeV. 1 

The general structure of the OPE expressions for moments is described elsewhere 
[15,21]. Here we give some technical details of our calculations. 
Charged lepton energy moments are computed as the ratios 

Mi n) (E cut ) = lEcut E 'F< . (2) 
If &dEe 

J-c-cut dEi 

Both numerator and denominator have been computed through order l/m\ compared to 
the leading partonic result; perturbative corrections are included from terms a s and floa 2 , 
and the ratio has been expanded in both perturbative and nonperturbative corrections. 
Since perturbative corrections in the Wilsonian scheme are suppressed, an alternative 
procedure with separate numerical evaluation of both numerator and denominator would 
yield a close result. Perturbative effects are expressed in terms of a^ is (m b ) for which we 
adopt 0.22 as central value. 

Hadronic invariant mass squared M\ in the OPE appears as a special choice in a 
one-parameter family of kinematic variables 

M 2 (L U ) = (L y +m b v u -q u ) 2 = L 2 + {m b v u -q u ) 2 + 2L v {m b v v -q v ) , 

{m b v v -q u ) 2 = m 2 x , (m b -q ) = e x (3) 



1 Here we use, however, the 'pole-type' Darwin expectation value p z D ~ p%(l GeV) — 0.1 GeV 3 instead 
of the Wilsonian p 3 D {l GeV), cf. Ref. [18]. 
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corresponding to L u = A'v u , where v u = ^ = (l,0) is the B four-velocity, and m 2 and e x 
have the meaning of invariant mass squared and energy in the hadron sector at the quark 
decay level, e x = Ex— A'. The OPE computes m 2 and e x much in the same way as lepton 
energy moments; therefore averages of M\ and its powers are generally polynomials in 
(M B —m b ), with the coefficients computed in the local OPE [22]. The latter are various 
mixed moments {e x m 2n ) derived from Eq. (1). 

The variable M 2 represents an observable for arbitrary L v . Being a combination of 
hadronic invariant mass and hadronic energy (and most generally of spacelike momentum), 
it can also be viewed as conventional "hadronic" invariant mass square if one considers 
not the decay of an isolated B meson, but rather of a compound of B combined with a 
non- interacting 'spurion' particle having momentum L v — {P^—mi ) v u ). From the OPE 
perspective, however, the inclusive probabilities appear as in the decay of the b quark, 
while the momentum associated with the light cloud looks like a spurion. The native 
object for the OPE is therefore M 2 (L U ) with vanishing L v , rather than L u = A'v u which 
yields M\. 

It turns out that for higher hadronic moments the generalized moments with L^-cA' 
are advantageous, they are better controlled theoretically and more directly sensitive to 
higher-dimension expectation values. 2 To utilize this advantage we consider, along with 
M\ the modified moments, i.e. those of 

M\ = M\- 2AE X + A 2 (4) 

with A = 0.65 GeV, close to the anticipated value of A'. In this case M\—m\ approximates 
the quark virtuality, for which higher moments with respect to average are intrinsically 
related to higher- dimension expectation values. (The last constant term in Eq. (4) does 
not affect such moments.) In our approach computing these modified moments does not 
require a new analysis - they are given by simply replacing M# by Mb — A. 

The higher moments of the decay distributions - in particular those of the hadronic 
invariant mass - are more informative when considered with respect to the average, say 
{{M x — (M x )) 2 ), {(M X —{M X )) 3 ), or similar moments for M x . These moments are the focus 
of our study. Since this may complicate to some extent the experimental error analysis, 3 
we also present similar numerical results for the moments evaluated with respect to a 
fixed hadronic mass, for which we take 4 GeV 2 in the case of M^-moments and 1.35 GeV 2 
for N"x ones. 

Power corrections in moments with lepton cut are obtain directly integrating the 
published heavy quark structure functions. The emerging analytic expressions are not 
too complicated, but lengthy and consist of many terms, especially for higher hadronic 
moments. We will present them in a dedicated publication [21]. Perturbative corrections 
are cumbersome and require numerical integrations. At the same time, in our approach 

2 This roots to lower infrared sensitivity of the modified moments compared to those of Ex which, 
in higher orders, arc dominated by maximal-g 2 kinematics where the c quark is nearly at rest. The 
combinatorial factors for the terms A'e x in the conventional M x moments work in the same direction, 
being additionally enhanced by large value of A'. 

3 We thank O. Buchmueller for discussing this point and alternative options. 
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utilizing a 'hard' separation between short- and long-distance effects they are numerically 
quite small. (For instance, it follows from the tables in the Appendix that including 
perturbative corrections in (M x ) has the same impact as decreasing m& by 15MeV, or as 
decreasing /4 by 0.02 GeV 2 in ((M 2 X -(M 2 X )) 2 ).) Therefore we evaluate hadronic moments 
to the first order in a s using the value a s = 0.3 (average gluon virtuality in B decays is 
lower than m^). Moreover, in the present paper we evaluated perturbative shifts in the 
hadronic moments neglecting the cut on the lepton energy. This seems to be a legitimate 
approximation since at E cut = they are small and do not exceed the expected accuracy, 
which is limited by other neglected effects. This element will be improved in Ref. [21]. 

Unlike the case of lepton moments, for power-suppressed terms we do not use the 
expanded form of the ratios that form the moments, nor drop any power-like terms gen- 
erated by p 2 , fj,Q, pp and p 3 LS wherever they appear. In particular, this implies that we 
would not include Mfi-mj (or its analogue M# — nib — A for A/j^-moments) into counting 
powers of Aq C d- This is natural since that entry is external to the OPE for inclusive 
probabilities and it can take arbitrary values, both much larger and much smaller than 
Aqcd, as illustrated by the modified hadronic moments. Similarly, in the perturbative 
corrections terms like ~ ct s /U 2 , a s PG e ^ c - are n °t retained, but those ~ ct s A k without 
nonperturbative expectation values are legitimately kept for arbitrary power k. A 

In this note we present our numerical results in a simplified form, using the following 
reference values of the parameters 



and providing the coefficients for the linear extrapolation in these values from this base 
point. Such linearized extrapolations appear sufficiently accurate for reasonable values of 
the parameters. The tables summarizing our results are given in the Appendix. Yet this 
is no more than the simplest compact way of communicating our results (for instance, 
we drew plots using the complete expressions rather than interpolations). The numeri- 
cal evaluations for other central values, more accurate interpolating tables, or compact 
Mathematica or FORTRAN programs evaluating them are available upon request, and 
will be provided with Ref. [21]. 

Measured in experiment are also non-integer moments of M x , most notably (M x ) and 
(M|) [23]. They do not arise naturally in the l/m& expansion, as illustrated by the limit 
m c ^0 which is analogous to the decay B — > X s + 7; fractional photon energy moments are 
not given there by the expectation values of local heavy quark operators. For B^X c £u, 
the OPE would involve an expansion in l/m c , as can also be seen from 



for integer moments with v = 2n the sum contains only terms through k — n and stops 
before (M|) enters the denominator. 

4 Terms 0(a s A' k ) with k>0 have not been included in the third hadronic moment. 



m b = 4.6GeV, m c =1.2GeV, p 2 =0.4GeV 2 , p 2 G = 0.35 GeV 2 , 
^ = 0.1 GeV 3 , pis = -0.15 GeV 3 , 



(5) 




(6) 
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Having computed three first (integer) hadronic moments, we truncated the sum in 
Eq. (6) after k = 3. Although incomplete, at the actual value of the charm mass this 
truncated expansion appears a sufficiently good numerical approximation for v — 1 and 
v = 3; the omitted terms seem significantly below the actual theoretical accuracy in 
evaluating the integer moments involved. 

In fact, M\ contains the dominant term 2A / e x as well as the perturbative bremsstrah- 
lung contribution, and both effects can be computed explicitly for arbitrary v without 
truncating the series in Eq. (6). The numerical impact of this resummation turns out 
totally negligible numerically. Moreover, while the term lh!e x dominates (M|-), it is 
subdominant in the higher moments with respect to (M|-}. 




Figure 1: Average hadronic invariant mass Figure 2: Second invariant mass moment 

squared (JW|) at different lepton energy ([M|- — (M|-)] 2 } (red) and second modified 

cuts, for the heavy quark parameters of hadronic moment {[N"x~ (-^"x)] 2 ) (blue), in 

Eq. (5). the same setting. 



For illustration we give a few plots showing dependence of the hadronic mass moments 
on the lepton energy cut, at the central values of heavy quark parameters. Fig. 1 depicts 
(M|), Fig. 2 addresses ((M| - (M|)) 2 ) and ((A/j| - (A/1)) 2 ). Non-integer M\ moments 
are shown in Fig. 3 where we actually plot the corresponding powers of the moments 
having dimension of mass, for v — 1, 2, 3, 4, 5 and 6. It shows that the full moments are 
by far dominated by the average invariant mass, with relatively small differences. 

As the semileptonic decay rate is well measured only above certain energy Ee, the 
OPE predictions for the decay fraction 

f t? ~irr dEp 

R*(E cut ) = JE p\* E * (7) 

are helpful to reconstruct the overall h^civ rate. Since the fraction of events cut out 
is small and they belong to the domain of low Ee which is theoretically most robust, the 
predictions for R* are expected to be quite reliable. 
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(M" x )i 




lepton energy cuts; heavy quark parameters as in Eq. (5). 



3 Discussion 



In the present paper we provide numerical expressions for the local OPE predictions of the 
charged lepton energy and hadronic mass moments with a lower cut on lepton energy. Two 
major aspects of our analysis have already been emphasized. First, no l/m c expansion is 
involved at any stage, and the c quark mass can more or less be arbitrary, large or small 
within the same formalism. The second aspect is that we rely on well-defined running 
quark masses and normalized at the scale 1 GeV. (The scheme is often referred to 
as 'kinetic', however this name is elucidating only when applied to quark masses.) The 
advantages of this approach are well known, and can be readily seen from the applications 
presented here. 

A quick comparison [24] with recent data [25-27] from DELPHI, BaBar and CLEO 
suggests agreement between this implementation of the OPE and experiment, with pre- 
ferred values of the heavy quark parameters in the range expected theoretically. (Our 
estimates show that the expectations for two moments of the photon energy spectrum 
appear to fit the values reported by CLEO [28] once the 'exponential' effects discussed 
in Ref. [24] are adjusted for.) The dependence of the first and second hadronic mass 
squared moments on the lepton energy cut seems to be in a qualitative agreement with 
the preliminary data reported by BaBar and CLEO [26,27]. The dedicated data analysis 
including possible fits to the predictions should, in our opinion be left for experiment, and 
we did not attempt this. 

Once the heavy quark parameters are extracted from data, one can readily determine, 
for example, \Va\ from the measured semileptonic b— >c decay width [18]: 

?S = (^lof 1 ) 2 (^) 2 (l-4.8[Br(S-X u M-0.0018])-[l + 0.30(a s (m fc )-0.22)] 

x [1 - 0.66 (m 6 (l GeV) -4.6 GeV) + 0.39 (m c (l GeV) - 1.15 GeV) 
+ 0.013 (/4 -0.4 GeV 2 ) + 0.09 (p 3 D (l GeV) -0.2 GeV 3 ) 

+ 0.05 (/4-0.35 GeV 2 ) -0.01 (p^+0.15 GeV 3 )] . (8) 
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It should be recalled, however, that p%(l GeV) appearing in the above equation is related 
to p 3 D employed in the present study by p 3 D ~ GeV)— 0.1 GeV 3 . 

We do not intend to address here the question of theory accuracy for the moments 
to an extent commensurate to the analysis [18] of total T s \(b — > c). Nevertheless, the 
numerical results presented in the Appendix allow a straightforward, if simplified, estimate 
of possible theory inaccuracies. 

Since no perturbative corrections to the Wilson coefficients of nonperturbative oper- 
ators have been computed so far, one can assume a related ~ 20% uncertainty in the 
contributions due to p^ and p? G and a ~ 30% uncertainty in those due to p 3 D and p 3 LS . 
As perturbative corrections in hadronic moments are presently implemented only at first 
order in a s , the associated uncertainty can be assessed by varying a s in a reasonable 
range; since the actual short- distance corrections selected by our Wilsonian treatment are 
moderate, we may conservatively vary the effective a s between 0.2 and 0.45. To safeguard 
against possible accidental cancellations in perturbative corrections, one may assume an 
additional minimal uncertainty in mj iC of about 20MeV. A similar 'minimal uncertainty' 
in pi is probably around 0.02 GeV 2 . 

The results we provide in the Appendix can be improved in several ways. The rele- 
vance of the improvements is essentially determined by the state of experiment. A more 
complete calculation of perturbative corrections in hadronic moments fully incorporating 
cuts in the lepton energy is the first on the list. All-order BLM resummation can be 
implemented in both lepton and hadron moments. There are ways to partially incorpo- 
rate non-BLM second-order corrections without performing extensive new calculations. 
All these improvements in the perturbative corrections are not expected to essentially 
modify the final result in the Wilsonian scheme (see, e.g. the dedicated analysis of T s \ in 
Ref. [18]). Yet having them incorporated would add confidence in the present estimates 
and may allow to reduce the theoretical uncertainty. 

A probe of the significance of second- and higher-order non-BLM corrections is already 
available by varying the normalization scale p used for quark masses and p^, p 3 D , while 
simultaneously running their values according to [29] 

K + !)-3(f-I)]) 

f (l«| + l)-3(f-i)]). 

H ln f + D- 3 Of -3)]). 

Hkg + iHOf-s)]), < 9 > 

where M denotes the scale used for normalizing a s . In fact, in the future we will also 
explicitly adopt the Wilsonian p 3 D (p) instead of p 3 D currently used. Even the uncertainty 
obtained in this way may, however, not fully represent the uncalculated non-BLM cor- 
rections. A further option is using the MS scheme for the charm mass, m c {m c ). With 
this choice, one could vary p in a wider range to probe more quantitatively the actual 
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hardness of a particular moment at a given Eg cut, and a more direct comparison with m c 
extracted from different physical processes [30] would be possible. Since no constraints is 
imposed on m b — m c through conventional mass relations like 

mfe - mc = — i i — + M^r *J + ° WJ + - (10) 

the normalization point - or even the scheme itself - can be different for m c and m b . 

It should be emphasized that we have quoted in the Appendix what literally emerges 
from our OPE expressions, regardless of whether the cut on Eg is mild or severe. However, 
if the cut is high, the effective hardness of the inclusive process degrades and the accuracy 
deteriorates. As has been argued [20,24], the truncated expressions of practical OPE may 
not fully reflect this, while the so-called 'exponential' effects may become significant when 
E cut exceeds f .3 to f .5 GeV. The same reservation applies to the above simplified way to 
estimate theoretical accuracy - it largely leaves out this aspect. The additional cut-related 
theory errors are possibly insignificant below 1.35 GeV; they may not dominate even at 
-^cut = 1-5 GeV - yet this cannot be confidently derived from theory a priori. The safest 
way to tackle them is to use as low E cut as possible, not larger than 1.2 GeV. 

A potential source of additional nonperturbative corrections are the so-called 'intrinsic 
charm', or IC contributions in the OPE, associated with non- vanishing expectation values 
of four-quark operators with charm fields bTc cTb. No allowance is presently reserved for 
them. It has been noted in Ref. [18] that their effect can mimic to some extent the 
contribution of the Darwin operator. We think therefore that at the moment theoretical 
constraints on p 3 D should be used cautiously in the context of fits based on the OPE 
without possible IC contributions. 

There are some conclusions we draw from our analysis. The moments under consid- 
eration seem to reliably (over) constrain a certain combination of the two heavy quark 
masses, approximately m^— 0.6m c ; the sensitivity to the individual masses is not high and 
they are subject to larger theoretical uncertainties. 

The moments are only weakly sensitive to the spin-orbit expectation value p 3 LS , so it 
cannot be extracted from the data. On the other hand, p 3 LS is reasonably constrained by 
a number of exact heavy quark sum rules which place model-independent bounds. Once 
they are taken into account, the associated uncertainty in the moments appears by far 
subdominant. Therefore, p 3 LS does not have to be included in the fit and we suggest to 
use a fixed value p\ s = —0.15 GeV 3 , and to vary it within ±0.1 GeV 3 to conservatively 
estimate the related uncertainty. This assumption would be further reinforced if the fits 
of data prefer values of the primary heavy quark parameters in the expected ranges 

m b {\ GeV) = 4.60 ± 0.06 GeV, /4(1 GeV) = 0.45 ± 0.1 GeV 2 , 

< p 3 D < 0.15 GeV 3 or 0.1 GeV 3 < p 3 D (l GeV) < 0.25 GeV 3 , (11) 

see Refs. [4,17,31]; an updated review of different determinations of beauty and charm 
quark masses can be found in Ref. [30]. Similarly, since p? G is accurately known, it should 
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be set to 0.35 GeV 2 [32] and varied within ±20% to allow for the perturbative uncertainty 
in its Wilson coefficient. 

Obtaining informative model-independent constraints on higher- dimension heavy quark 
parameters requires measurement of higher hadronic moments, at least the second and 
desirably the third one. Since these higher moments - when considered with respect to the 
average - depend strongly on higher- dimension expectation values, even a rough measure- 
ment of the variance and of the asymmetry parameter would yield precious information. 
One should be warned, however, that the theoretical accuracy one can realistically achieve 
for the higher moments of M x is limited, especially in the case of the third moment. On 
the other hand, the modified higher hadronic moments ((A/j| — (A/"|-)) 2 ), (A/j^-)) 3 ) 
are better in this respect and more suitable to this purpose. We encourage experiments 
to pay closer attention to such combinations of M\ and Ex moments. 

The actual assessment of the theoretical error in the calculation of the moments is a 
very subtle issue. We cannot do without mentioning a few important aspects. 

First, there are strong correlations among theoretical uncertainties. After all, every- 
thing inclusive we compute is expressed in terms of only three structure functions which 
have physical properties like positivity, regardless of any dynamics. Because of these cor- 
relations, one should distinguish between the overall consistency of the fits to the various 
moments in the context of our formalism, and the concrete prediction for a certain mo- 
ment. There are common uncertainties, like a systematic bias in or /j,q, that would 
simply shift the fitted value of m b or {1q, but would not degrade the quality of the global 
fit, nor would they alter significantly the -E cut -dependence [20]. At the same time, these 
uncertainties could still noticeably affect the numerical value of a particular moment. 

On the other hand, there are uncertainties that affect each moment in a different way, 
for instance, unknown perturbative corrections to the Wilson coefficients of the nonper- 
turbative operators, which are in principle different for every moment, and also depend on 
E cut . Therefore, simply varying the values of the heavy quark parameters uniformly in all 
observables may not represent realistically the uncertainty of the theoretical expressions, 
and an intermediate procedure may be required. 

Having mentioned these complications of the theory error analysis, we would like to 
make a few suggestions that can be inferred from our study. For sufficiently low cut on 
En a reasonable starting point is to estimate the theoretical accuracy in the moments just 
varying the values of the heavy quark parameters they depend upon in the ranges we have 
mentioned above Eq. (9). For higher moments, however, this should be applied to the 
moments with respect to average; the ordinary moments (around zero) would in this way 
exhibit strong theory error correlations. Clearly, allowance should be made for additional 
uncertainty once the cut is raised beyond 1.35 GeV. 

In summary, we have seen that consistency checks between theory and emerging data 
should rely on robust theoretical elements, with additional assumptions reduced to mini- 
mum. Once agreement with theory is confirmed, and the domain and degree of applica- 
bility is verified experimentally (e.g., the safe interval for lepton cut), one can and should 
implement into the fit of data all the model-independent relations and bounds follow- 
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ing from heavy quark sum rules. These are exact relations derived in QCD, and they 
are indispensable for a precision determination of the heavy quark parameters and for a 
stringent test of our theoretical tools. 
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Appendix 

The following tables give our numerical estimates for various moments. The general form 
for a generic moment M. in question will be 

M (m b , m c , nl , fia > Pd , Pls ; a a ) = V + B (m b - 4.6 GeV) + C (m c - 1 .2 GeV) 

+ P(^_o.4GeV 2 ) +D{p 3 D -0.1 GeV 3 ) (A.l) 
+ 0.35 GeV 2 ) + L (p^+0.15 GeV 3 ) + S(a s -0.22) ; 

V represent the reference values obtained for the heavy quark parameters in Eq. (5). 
They have dimension of the moment M. itself and the quoted number is in GeV to the 
corresponding power. The values of all the coefficients B to S are likewise in the proper 
power of GeV (the same power for S, one power less for B and C, two powers less for P 
and G, three powers less for D and L). Values of E cut are shown in GeV as well. 

All the moments are given without cut on E e (i.e. E cut = 0) and for E cut = 0.6 GeV, 
E cnt = 0.9 GeV, -E cut = 1-2 GeV and E cut = 1.5 GeV. Fraction ratios R* are tabulated for 
E cut = 0.3 GeV, E cut = 0.6 GeV and E cut = 0.9 GeV. 

A.l Lepton energy moments 



Table 1. First moment of the lepton energy (Eg). 





V 


B 


C 


P 


D 


G 


L 


S 





1.372 


0.389 


-0.286 


0.033 


-0.085 


-0.078 


0.0043 


-0.029 


0.6 


1.420 


0.367 


-0.268 


0.036 


-0.083 


-0.076 


0.0035 


-0.032 


0.9 


1.500 


0.339 


-0.243 


0.041 


-0.081 


-0.073 


0.0021 


-0.035 


1.2 


1.614 


0.304 


-0.211 


0.051 


-0.083 


-0.068 


-0.0005 


-0.043 


1.5 


1.759 


0.262 


-0.173 


0.073 


-0.096 


-0.063 


-0.0058 


-0.055 
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Table 2. Second moment of the lepton energy with respect to average, {(Eg — (Eg)) 2 } 



E C ut 


V 


B 


C 


P 


D 


G 


L 


S 





0.1774 


0.0951 


-0.0584 


0.044 


-0.059 


-0.029 


-0.0048 


-0.036 


0.6 


0.1391 


0.0960 


-0.0602 


0.042 


-0.058 


-0.027 


-0.0047 


-0.033 


0.9 


0.0969 


0.0868 


-0.0547 


0.040 


-0.057 


-0.024 


-0.0048 


-0.031 


1.2 


0.0570 


0.0691 


-0.0433 


0.038 


-0.055 


-0.019 


-0.0051 


-0.029 


1.5 


0.0266 


0.0473 


-0.0293 


0.035 


-0.057 


-0.012 


-0.0056 


-0.028 



Table 3. Third moment of the lepton energy with respect to average, ((Eg — (Eg)) 3 ) 





V 


B 


C 


P 


D 


G 


L 


100-5 





-0.0334 


-0.0307 


0.0265 


0.023 


-0.034 


0.0043 


-0.0054 


-1.6 


0.6 


-0.0121 


-0.0204 


0.0189 


0.022 


-0.034 


0.0021 


-0.0049 


-1.7 


0.9 


-0.0015 


-0.0082 


0.0086 


0.019 


-0.033 


0.0003 


-0.0041 


-1.5 


1.2 


0.0016 


0.0001 


0.0015 


0.015 


-0.030 


-0.0005 


-0.0031 


-1.2 


1.5 


0.0009 


0.0030 


-0.0013 


0.010 


-0.025 


-0.0003 


-0.0020 


-0.9 



Table 4. Fraction of the decay rate R* with Eg exceeding a threshold value E cut 



-^cut 


V 


B 


C 


P 


D 


G 


L 


100-5 


0.3 


0.9934 


0.006 


-0.0044 


-0.0001 


-0.0008 


-0.0008 


0.0001 


0.06 


0.6 


0.9508 


0.040 


-0.0317 


-0.0008 


-0.0060 


-0.0060 


0.0009 


0.15 


0.9 


0.8476 


0.120 


-0.0947 


-0.0018 


-0.0176 


-0.0189 


0.0029 


0.19 



A. 2 Hadron invariant mass moments 

Since the moments related to hadronic invariant mass and hadronic energy are presently 
calculated using only 0(a s ) perturbative corrections, we have employed a s = 0.3, which 
can equivalently be represented as 

»s(m b ) (A2) 



1_ 0.846|^ 

2 7T 

with the canonical value a s (m^} = 0.22 . This allows to use the same form Eq. (A.l) 
with S showing the sensitivity to perturbative corrections. In this way the interval 0.2 < 
a s <0.45 mentioned in Sect. 2 corresponds to varying a s (m b ) within 0.22l;2 6 . Switching 
off perturbative corrections, a s = 0, numerically amounts to subtracting 0.16 5 from a 
moment. 

As seen from the tables, the dependence on the precise value of a s is moderate and the 
corresponding uncertainty is subdominant; this is an advantage of using the Wilsonian 
scheme. 
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Table 5. First hadronic invariant mass moment (M^-). S'~0.41 



^cut 


V 


B 


C 


P 


D 


G 


L 





AMI 


-4.99 


3.18 


-0.70 


1.0 


0.48 


-0.13 


0.6 


4.619 


-4.96 


3.18 


-0.73 


1.0 


0.50 


-0.13 


0.9 


4.583 


-4.90 


3.19 


-0.79 


1.1 


0.55 


-0.12 


1.2 


4.528 


-4.84 


3.21 


-0.93 


1.3 


0.64 


-0.11 


1.5 


4.444 


-4.74 


3.22 


-1.29 


1.8 


0.81 


-0.06 



Table 6. Second invariant mass moment with respect to average, ((M|-(M|)) 2 ). S 1 ^— 0.51 



^cut 


V 


B 


C 


P 


D 


G 


L 





1.290 


0.396 


-0.97 


4.8 


-5.9 


-0.14 


0.30 


0.6 


1.233 


0.446 


-0.98 


4.7 


-6.0 


-0.13 


0.24 


0.9 


1.137 


0.523 


-0.99 


4.5 


-6.4 


-0.11 


0.13 


1.2 


0.985 


0.660 


-1.00 


4.2 


-7.1 


-0.06 


-0.05 


1.5 


0.747 


0.970 


-1.04 


3.8 


-8.9 


0.05 


-0.34 



Table 7. Third invariant mass moment with respect to average, ((Mj| — (M^)) 3 }. 5 ~ 10.1 



^cut 


V 


B 


C 


P 


D 


G 


L 





5.02 


1.40 


-2.50 


5.6 


21 


-1.3 


1.2 


0.6 


4.92 


1.48 


-2.49 


5.5 


21 


-1.3 


1.1 


0.9 


4.80 


1.53 


-2.45 


5.5 


19 


-1.3 


1.0 


1.2 


4.68 


1.52 


-2.32 


5.7 


17 


-1.2 


0.8 


1.5 


4.64 


1.40 


-1.98 


6.3 


15 


-0.8 


0.4 



Table 8. First modified hadronic moment (7VJ-). ,S~0.38 



^cut 


V" 




C 


P 


D 


G 


L 





1.580 


-4.06 


2.52 


-0.62 


0.74 


0.34 


-0.10 


0.6 


1.579 


-4.03 


2.53 


-0.64 


0.77 


0.36 


-0.10 


0.9 


1.574 


-3.98 


2.53 


-0.69 


0.83 


0.40 


-0.09 


1.2 


1.558 


-3.90 


2.52 


-0.82 


0.96 


0.47 


-0.08 


1.5 


1.501 


-3.75 


2.45 


-1.15 


1.30 


0.57 


-0.04 
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Table 9. Second modified hadronic moment with respect to average ((A/j|— (A/j|)) 2 ). S~— 0.35 



E C ut 


V 


B 


C 


P 


D 


G 


L 





0.954 


0.506 


-0.606 


3.6 


-3.8 


0.12 


0.25 


0.6 


0.918 


0.537 


-0.616 


3.5 


-3.9 


0.12 


0.20 


0.9 


0.855 


0.592 


-0.624 


3.4 


-4.1 


0.14 


0.13 


1.2 


0.761 


0.685 


-0.628 


3.2 


-4.6 


0.18 


0.00 


1.5 


0.621 


0.882 


-0.659 


3.0 


-5.8 


0.31 


-0.17 



Table 10. Third modified hadronic moment with respect to average {{M x — (A/j|)) 3 ). S*~9.9 



E cu t 




B 


C 


P 


D 


G 


L 





3.27 


1.47 


-1.67 


2.5 


13.3 


-0.87 


0.32 


0.6 


3.22 


1.53 


-1.68 


2.4 


12.8 


-0.88 


0.29 


0.9 


3.15 


1.57 


-1.66 


2.5 


12.1 


-0.89 


0.25 


1.2 


3.10 


1.57 


-1.58 


2.7 


11.1 


-0.88 


0.17 


1.5 


3.16 


1.43 


-1.31 


3.6 


9.5 


-0.69 


-0.03 



Tables 11 and 12 give predictions for non-integer moments (Mx) and (M x ); they are 
evaluated using Eq. (6) truncated after k = 3. 



Table 11. Hadronic invariant mass moment (Mx)- S~0.12 



^cut 


V 


B 


C 


P 


D 


G 


L 





2.145 


-1.17 


0.750 


-0.22 


0.33 


0.11 


-0.03 


0.6 


2.140 


-1.17 


0.752 


-0.22 


0.34 


0.12 


-0.03 


0.9 


2.133 


-1.16 


0.757 


-0.24 


0.37 


0.13 


-0.03 


1.2 


2.122 


-1.15 


0.765 


-0.27 


0.42 


0.15 


-0.02 


1.5 


2.105 


-1.14 


0.770 


-0.35 


0.56 


0.19 


-0.01 



Table 12. Hadronic invariant mass moment (M x ). 



^cut 


V 


B 


C 


p 


D 


G 


L 





10.19 


-16.0 


10.1 


-1.4 


2.1 


1.5 


-0.38 


0.6 


10.11 


-15.8 


10.1 


-1.5 


2.1 


1.6 


-0.38 


0.9 


9.98 


-15.6 


10.1 


-1.8 


2.3 


1.7 


-0.38 


1.2 


9.78 


-15.2 


10.1 


-2.2 


2.8 


2.0 


-0.37 


1.5 


9.47 


-14.8 


10.0 


-3.4 


3.9 


2.6 


-0.27 



Below we give the tables for the higher integer hadronic moments (for mass squared 
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and modified) with respect to a fixed mass. 



Table 13. Second invariant mass moment with respect to 4GeV 2 , ((M^ — 4 GeV 2 ) 2 ) . 5 



E C ut 


V 


B 


C 


P 


D 


G 


L 





1.701 


-5.76 


3.20 


3.9 


-4.6 


0.47 


0.13 


0.6 


1.617 


-5.45 


3.06 


3.8 


-4.7 


0.50 


0.08 


0.9 


1.477 


-4.95 


2.83 


3.6 


-5.0 


0.54 


-0.02 


1.2 


1.264 


-4.21 


2.49 


3.2 


-5.7 


0.62 


-0.17 


1.5 


0.944 


-3.02 


1.92 


2.7 


-7.4 


0.78 


-0.40 



Table 14. Third invariant mass moment with respect to 4 GeV 2 , ((M| — 4 GeV 2 ) 3 ). S~ll 





V 


B 


C 


P 


D 


G 


L 





7.77 


-22.9 


12.0 


11.2 


15.0 


0.85 


1.09 


0.6 


7.45 


-21.3 


11.2 


10.6 


14.2 


0.88 


0.91 


0.9 


6.98 


-18.9 


10.0 


9.7 


12.8 


0.96 


0.63 


1.2 


6.39 


-15.5 


8.3 


8.7 


10.8 


1.15 


0.25 


1.5 


5.72 


-10.6 


5.8 


7.6 


7.5 


1.62 


-0.27 



Table 15. Second modified moment with respect to 1.35 GeV 2 ((A/1-1.35 GeV 2 ) 2 ). 5 ~ -0.19 



^cut 


V 


B 


C 


P 


D 


G 


L 





1.007 


-1.20 


0.62 


3.3 


-3.4 


0.28 


0.21 


0.6 


0.970 


-1.14 


0.60 


3.2 


-3.5 


0.29 


0.16 


0.9 


0.905 


-1.03 


0.57 


3.1 


-3.7 


0.32 


0.08 


1.2 


0.805 


-0.78 


0.49 


2.9 


-4.2 


0.38 


-0.03 


1.5 


0.644 


-0.11 


0.14 


2.7 


-5.4 


0.49 


-0.19 



Table 16. Third modified moment with respect to 1.35 GeV 2 , ((M% - 1.35 GeV 2 ) 3 ). 5~11 



^cut 


V 


B 


C 


P 


D 


G 


L 





3.94 


-10.4 


5.5 


3.0 


12.8 


0.25 


0.19 


0.6 


3.86 


-9.8 


5.2 


2.9 


12.3 


0.26 


0.15 


0.9 


3.74 


-8.8 


4.8 


2.8 


11.5 


0.28 


0.08 


1.2 


3.59 


-7.4 


4.1 


2.7 


10.4 


0.36 


-0.02 


1.5 


3.45 


-5.5 


3.1 


2.6 


9.1 


0.56 


-0.18 



5 There are strong cancellations in the perturbative corrections here, and the face value of S is not 
meaningful. One can roughly assume |S'|^0.25. 
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